Current Induced Order Parameter Dynamics: Microscopic Theory Applied to 

Co/Cu/Co spin valves 
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Transport currents can alter alter order parameter dynamics and change steady states in super- 
conductors, in ferromagnets, and in hybrid systems. In this article we present a scheme for fully 
microscopic evaluation of order parameter dynamics that is intended for application to nanoscale 
systems. The approach relies on time-dependent mean-field-theory, on an adiabatic approximation, 
and on the use of non-equilibrium Greens function (NEGF) theory to calculate the influence of a 
bias voltage across a system on its steady-state density matrix. We apply this scheme to examine the 
spin-transfer torques which drive magnetization dynamics in Co/Cu/Co spin- valve structures. Our 
microscopic torques are peaked near Co/Cu interfaces, in agreement with most previous pictures, 
but suprisingly act mainly on Co transition metal d-orbitals rather than on s-orbitals as generally 
supposed. 
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PACS numbers: 71.20.Lp,71.20.Be,61.50.Lt 



I. INTRODUCTION 



Transport currents can be used to alter the order pa- 
rameter dynamics of metals with broken symmetries, 
including ferromagnets, antiferromagnets, superconduc- 
tors, and hybrid systems containing both magnetic and 
superconducting elements. A familiar and simple exam- 
ple of this type of phenomena is a circuit in which current 
flows between normal (N) metal elements though a su- 
perconductor (S). Because of Andreev reflection of the 
quasiparticle current at the N/S boundary, the Cooper 
pair amplitude in the superconductor is altered. In the 
presence of a normal metal transport current, the super- 
conducting order parameter developes a spatial gradient 
in the steady state which allows the condensate to carry 
current through the superconducting element. A more re- 
cent example is quantum Hall bilayers which develop ex- 
citonic condensates^ under certain conditions and exhibit 
anomalies related^ to the influence of transport currents 
on interlayer phase coherence. An important class of re- 
lated phenomena which has received a large amount of 
theoretical^iii^i^ and experimenta&^i^ attention over the 
past decade is centered on the influence of transport cur- 
rents on magnetization in magnetic metals, particularly 
ferromagnetic nanoparticles and ferromagnets containing 
domain walls. In this case the current-related torques 
exerted on the magnetization can usually be understood 
at least qualitatively by appealing to conservation of to- 
tal spin: the magnetization torques arc understood as 
the reaction counterpart to the action of the magnet's 
collective exchange field on the quasiparticles, i.e. as 
spin-transfer torques (STTs). In STT physics, transport 
electrons can change the magnetic state of the device, 
by switching the orientation, or inducing microwave fre- 



quency oscillations in the orientation of magnetic layers. 
These phenomona hold out the promise of applications, 
for example for writing magnetically stored information. 
Calculations of STTs generally proceed by computing the 
spatial dependence of spin currents in a circuit and invok- 
ing conservation of spin angular momentum to infer the 
torque which acts on the magnetization. This approach 
has so far had mixed levels of success in predicting im- 
portant quantities like critical currents for magnetization 
switching. Most calculations to date have been based on 
free electron or on semi-empirical tight-binding models. 

In this paper we discuss a practical scheme for es- 
timating the influence of transport currents on order- 
paramater dynamics in superconductors, magnetic met- 
als, or hybrid systems. Our approach, which is ex- 
plained in detail in Sec. [Ill starts by assuming a time- 
dependent mean-fleld theory in which the electron dy- 
namics is described by a single-particle Hamiltonian 
which is uniquely deflned by the density-matrix of the 
system. In practice the most flexible and powerful choice 
for such a time-dependent mean-field theory will nor- 
mally be spin-density functional theory, although any 
similar mean-field-like approximation is consistent with 
the procedures outlined below. Hartree-Fock theory, for 
example, has the same structure and has the advan- 
tage of being able to describe non-local exchange effects 
which might be important under some circumstances, 
for example when current is carried by a partially occu- 
pied bonding band. We then make an adiabatic approx- 
imation by assumming that the time-dependence of the 
single-particle Hamiltonian can be ignored, and use non- 
equilibrium Greens function (NEGF) theory to evaluate 
the influence of a bias voltage on the system density ma- 
trix. Our implementation of NEGF theory is explained in 



Sec, mil The adiabatic approximation is normally appro- 
priate since the driving terms are proportional to trans- 
port current and are typically weak compared to charac- 
teristic energy scales. Since our interest is in nanoscale 
systems, the neglect of inelastic scattering in the system 
implied by the use of steady state NEGF theory is nor- 
mally justified. (Inelastic scattering lengths under ambi- 
ent conditions are typically ~ lOnm. The approach could 
be extendediS, to include inelastic scattering.) The influ- 
ence of current on order parameter dynamics follows from 
the dependence of the density-matrix on bias voltage for 
any given single-particle Hamiltonian. 

In the remaining sections of the paper, we apply our 
theory of current-induced order parameter dynamics to 
a Co/Cu/Co spin- valve system, an example of transition 
metal ferromagnet spintronics which is of great practical 
importance and has recieved considerable attention. In 
Sec. IIVI we descibe this ab initio calculation, for which 
we choose a geometry similar to that studied in Ref. [l^, 
in detail. This calculation is intended to illustrate that 
ab initio detail will play a key role in designing mag- 
netic nanostructures whose non-linear magnetotransport 
properties are optimized for particular applications. In 
Sec. lYl we review the results, and discuss some of the 
surprising insights that this approach gives into the mi- 
croscopic physics of spin transfer. Finally, in Sec. IVI) 
we review our results and outline some other interesting 
potential applications of this microscopic theory of order 
parameter dynamics in nanostructures. 



II. MICROSCOPIC THEORY OF 
CURRENT-INDUCED ORDER PARAMETER 
DYNAMICS 



The practical scheme for microscopic calculation of 
order parameter dynamics outlined in this section is 
based on a previously decribed microscopic theory of 
spin-transfer torques in circuits containing ferromagnetic 
metal nanoparticles?ii This scheme views spin-transfer 
torques as a specific example of a more general class of 
phenomena in which collective dynamics is altered by 
a bias voltage because of the change it induces in the 
relationship between the single-particle density matrix 
and the effective single-particle Hamiltonian. This micro- 
scopic view of current-induced order parameter dynamics 
in magnetic metals has consequences that might be signif- 
icant in some instances, for example in circuits contain- 
ing antiferromagnetic metal elements?^ In this section 
we explain the approach, using a notation that is con- 
venient for the NEGF calculations we apply Co/Cu/Co 
spin-valves in the body of the paper. 



A. Order Parameter Dynamics in an Isolated 
System 

When an interacting electron system is described by 
a time-dependent mean-field theory the effective single- 
particle Hamiltonian can be constructed at each instant 
in time from the single-particle density matrix. His- 
torically the first example of this type of theory is the 
Hartree-Fock approximation. In modern usage density- 
functional-theory is usually more accurate and easier 
to implement. For example, in the local spin-density- 
approximation, the Hamiltonian can be constructed from 
the density and the spin-density, i.e. from elements of 
the single-particle density matrix p that are diagonal in 
a position representation. The exchange and correlation 
potentials in the Hamiltonian H are explicit non-linear 
functions of the density matrix: H = H[p]. 

We assume that a procedure for constructing the single 
particle Hamiltonian from the density matrix has been 
chosen. Then the dynamics of p is specified by its equa- 
tion of motion: 

a,p=l[H[p],p]. (1) 

For systems with broken symmetries, for example ferro- 
magnets, the order parameter is specified by a particular 
average of the density matrix so that Eq. ^ specifies its 
dynamics. In this purely electronic equation of motion we 
have ignored coupling between the electronic system and 
nuclear spins, phonons, or other environmental degrees 
of freedom which can play a role in some circumstances. 
Such couplings arc accounted for by phcnomcnological 
additions to the right hand side of Eq. ^ or to appro- 
priate averages of this equation. 

This paper focuses on the influence of transport cur- 
rents produced by a bias voltage on the collective mag- 
netization dynamics of a small magnetic nanoparticle. It 
is helpful to consider first an isolated system in which 
transport currents are absent. We separate both the 
single-particle Hamiltonian and the density matrix into 
its spin-dependent and spin-independent contributions: 

P^'^■s's - ^[P,,,ds>s+m,H-Ts>s\. 

where r is the vector of Pauli spin matrices, i'i are or- 
bital indices, and s's are spin-indices. The notation for 
the spin-dependent part of the Hamiltonian is chosen to 
emphasize that it produces a spin-splitting A when it 
is orbital independent, as often assumed in simple toy 
models of a ferromagnetic metal. In mean field approx- 
imations, the interaction contribution to A and m are 
locally related according to 

-* Tn 

A = Ao(n,m)- (3) 
m 

where n and m are the local charge and spin densi- 
ties, respectively, and Aq is some paramerization of the 
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exchange-correlation potential. Given this notation, it 
is possible to derive a useful expression for the time- 
dependence of the the a-th component of the spin-density 
5 of a chosen subsystem (SS) 
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(4) 



The first four terms on the right hand side of Eq. ^ rep- 
resent the net spin-current into the subsystem which has 
contributions from both the spin-polarization of inter- 
orbital coherence and from the spin-dependence of the 
inter-orbital matrix elements in the Hamiltonian. In 
the final term, contributions in which j ^ SS represent 
an additional spin current from spin-dependent hopping. 
The contributions from j G SS describes precessional 
time-evolution of spins in the subsystem under the in- 
fluence of effective magnetic fields implied by the spin- 
dependent terms in the Hamiltonian. We make use of 
this expression below. 

It follows from Eq. ([4]) that the total spin of an isolated 
magnetic nanoparticle satisfies 



= ^Tr[pf] = ^Tr[m x A] 



(5) 



If we assume that the magnetization in the nanomagnet is 
coUinear (or at least nearly so), the interaction contribu- 
tion to the spin-dependent part of the Hamiltonian, will 
be in the same direction as the magnetization, according 
to Eq. ([3]). It follows that electron-electron interactions, 
which always provide the dominant contribution to A, 
do not directly influence the total spin-dynamics. When 
the only other spin-dependent terms in the Hamiltonain 
are due to an external magnetic field which produces an 
orbital independent splitting field A^'^^ Eq. ^ reduces 
to 



hS = SxA'^^l 



(6) 



More generally, additional spin-dependent terms in the 
Hamiltonian due to spin-orbit coupling and magnetic 
dipole-dipole interactions produce additional contribu- 
tions to the effective field so that A''^' is replaced by 
an effective magnetic field A^'^-'''^) that depends on the 
total spin orientation and includes magnetocrystalline 
anisotropy, shape anisotropy, and the dissipative contri- 
bution due to the coupling between the total spin and in- 
coherent particle hole excitations. These effects are nor- 
mally described phenomenologically using the Landau- 
Liftshitz Gilbert equation: 



hS = Sx A^^ff^ + ah-^ X S. 

\s\ 



(7) 



B. Influence of a Bias Voltage 

We now address corrections to the Landau-Liftshitz- 
Gilbert equations which apply when a bias voltage is 
applied to a metallic nanomagnet. We first describe 
our strategy for a first principles description of current- 
induced order parameter dynamics using a more general 
language with wider applicability and then specialize to 
the magnetic spin-transfer torque case. Our approach is 
based on the non-equilibrium Greens function description 
of non-interacting fermions under the infiuence of a bias 
voltagei^ which we describe in more detail below. In this 
approach, the bias voltage is represented by placing the 
system in contact with particle reservoirs with chemical 
potentials — + eVe/S and pd — ep — cVb/^ in 
source and drain respectively. When the system Hamil- 
tonian and its coupling to source and drain electrodes 
is time- independent, electrons with energies inside of 
the transport window ^ld < E < (is solve a time- 
independent Schroedinger equation with incident-from- 
source scattering boundary conditions. The Schroedinger 
equation solution for the system plus reservoirs system 
is readily constructed using Greens function techniques. 
Electrons inside and outside the transport window be- 
have very differently. Given a Hamiltonian, the density 
matrix may be separated into a contribution from elec- 
trons inside and outside the transport window: 



Ptot[H] = /9con<i[H] + Pf,.[H] 



(8) 



where we define pcon(i[H] as the contribution to the to- 
tal density matrix from all states with energy below the 
transport window minimum /Lt£), and pt, [H] as the con- 
tribution from states in the transport window. (The no- 
tation pcond is intended to suggest a condensate density- 
matrix since we will ultimately take an average of the 
density matrix which highlights an observable associated 
with a broken symmetry.) To make progress we limit our 
attention to circumstances in which the order parame- 
ter and hence the effective Hamiltonian changes slowly in 
time and the relevant contributions to the density matrix 
are dominated by orbitals outside the transport window. 
Given these approximations, which are soundly based for 
magnetic systems at least, 



Plot = i[H,pco«rf](l -r -Q^ 

^ Ptot « »[H, Pconcl\ 



(9) 
(10) 



Where we have noted that the second term of Eq. Q is 
of order eV/ep, which is negligibly small for metallic sys- 
tems. The current-induced contribution to the density- 
matrix equation of motion is due to the difference be- 
tween li[ptot] and Hlpcond]- 



PCI 



condl : Pc 



(11) 



Making the appropriate average for the total spin of a 
nanoparticle we obtain the version of Eq. ([5]) which de- 
scribes the current-induced contribution to the dynamics 
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of the total nanoparticle spin: 



1, 



= \'Tr[rncond x At^]. 



(12) 



where Ajot and l^cond are the spin-dependent Hamitoni- 
ans for the system arising from ptot and Pcond, respec- 
tively, and Atr is the difference between the two. As de- 
scribed above, we consider situations in which the states 
in the transport window have relatively little contribu- 
tion to the order parameter (rricond ^ TTiitr), and focus 
on exchange contributions to A, which arc dominant, so 
that, according to Eq. ([3]) 



Atot = Ao(n, m) 



[mcond + mtr) 



Ktr = Ao(ri, m) 



(13) 



To evaluate the current-induced magnetization dynam- 
ics, we need to evaulate the right hand side of Eq. (fT^ . 

We assume that we have a circuit containing a nano- 
magnet with approximately coUinear magnetization. For 
electrons in the transport window we can apply Eq. Q 
to obtain for any orbital i in the system 



= 



— y 



(0) 



(0) 



p(°)A" 1 



UNP 



2 (Aconrf X rhtr) 



(14) 



In the above equation, p^^^ and rh refer to the trans- 
port contribution to the total density matrix. In writ- 
ing the second term as a simple cross product, we have 
also made a simplifying assumption that spin-dependent 
hopping terms are negligible, which is valid for the sys- 
tems studied here. As explained earlier, the first group 
of terms on the right-hand-side of Eq. (fT4|) represents 
is the difference between the spin-current into and out 
of the nanoparticle (NP), while the last term represents 
the precession of transport electron spin in the presence 
of the exchange field from the condensate. In order to 
change their spin-polarization as they move from source 
to drain through the nanomagnet, transport electrons 
must align their spins at an angle to the exchange field. 
The exchange field produced by the condensate electrons 
produces a torque on the current carrying quasiparticles. 
Comparing with Eq. (fT^ . we see that this torque is equal 
and opposite to the torque applied to the condensate by 
the transport electrons. 

In summary, to find the current induced torque on the 
magnetization, we evaluate Eq. (|12p . or, equivalently, 
the 2nd term of Eq. ([M]). Eq. ([14]) implies that this 
torque is equal to the net spin current into or out of 
the nanoparticle, which is the standard picture of spin 
transfer torque, and which relies on conservation of total 



spin angular momentum arguments. As emphasized in 
Ref. [Ill, our more general approach is applicable to sit- 
uations in which there is no conservation of spin, and is 
therefore imperative in systems with spin-orbit coupling, 
and in anti-ferromagnetic systems, for example. 

Finally, to gain further insight into the physics of STT, 
we evalute Eq. (fT2|) locally by choosing the subsystem to 
be a single atom, or single atomic orbital, as opposed 
to the entire nanoparticle. To obtain the total STT on 
the nanoparticle, we sum over these individual contribu- 
tions. Such a partition of STT is strictly speaking only 
valid in the case where hopping is spin-independent 
which is approximately the case for our calculations. We 
have explicitly checked that the difference in STT on the 
nanoparticle found is this manner compared to evaluat- 
ing Eq. globally is negligible. 



III. NON-EQUILIBRIUM GREENS FUNCTION 
CALCULATIONS 



To implement the ideas in Sec. ITTJ we use non- 
equilibrium Green's functions (NEGF), within a den- 
sity function theory framework. Briefly, the Hamiltonian 
of the device is calculated within standard parameter- 
izations of LSDA)^ generalized to include noncoUinear 
magnetization. The electronic states are then popu- 
lated according the the distribution functions of the left 
and right semi-infinite leads, determining a density ma- 
trix p[H]. The details of the calculation procedure have 
been given previously. ^'^ Such formalism has been used 
to successfully calculate transport properties of similar 
systems r^ii^ii^ and spin-dependent properties of molec- 
ular scale systems 

Since the system under consideration here is in the 
metallic regime, we work within linear response approx- 
imation for the bias dependence Vb- In this case, the 
expression for the current induced torque on atom i is 
given as 



Sci 



id(i,i3) X mtr(j,/3))Vs- (15) 



where the (3 is an orbital label for atom i, and refers 
to the transverse momentum of each propagating state. 
The transport contribution to the density matrix is 



Ptr 



G''lm{T.l)G'' . 



(16) 



Here E£ is the retarded self energy, which accounts for 
the presence of left semi-infinite left lead, and G""'" are 
the retarded (advanced) advanced Green's function for 
the device. The above Green's functions and self energy 
are evaluated at the fermi energy, and the term in- 
dicates that non-equilibrium electrons emanate from the 
left lead. We have checked explicitly that the transport 
density matrix satisfies Eq. 
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Within the Landauer formaUsm (and staying within 
Unear response), the current density for a bias Vb is: 

J-j^ji^Il J dhn^.^A^^)yB. (17) 

where the transmission coefBcient is given by: 

Tk,.a,a' = Tr[(Im(I]£)G'-Im(S^j)G"),,,,]. (18) 

The trace above refers to orbital and site. Since most 
experiments are done under fixed current bias, an exper- 
imentally relevant quantity is the total spin torque per 
current, given by the ratio of Eq. P3|) (summed over all 
atoms) to Eq. (fT7|) . In linear response, this is: 

S _ HE J2t,i3i^cond{-i.,P) X mtr(t,i3)) 

IV. CALCULATION DETAILS 
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FIG. 1: Conductivity versus relative layer orientation. 



The system studied consists of a semi-infinite Co lead, 
a Cu spacer layer with 9 atomic planes, a Co free layer 
with 15 atomic planes, and a semi-infinite Cu lead. All 
layers have translational symmetry in the transverse di- 
rection, and therefore represent the bulk realistically. 
Both Co and Cu are assumed to be in the fee phase, and 
we use a lattice constant of 3.54 A throughout. We use 
norm-conserving pseudopotentials^^ and an s,p, d single- 
zeta basis set. We have found excellent agreement with 
established band structure, density of states, and bulk 
conductivity for Co and Cu with this basis set. We 
have found that 800 fc-points within the Brillouin zone 
is sufficient for convergence of the self-consistent density 
matrix. To calculate transmission coefficients and non- 
equilibrium spin densities, we have used 25,600 fc-points. 
We have found that there is less than a 1 % difference in 
these quantities when using up to 32,400 fc-points. In cal- 
culating the spin torque contribution for different chan- 
nels in the Brillouin zone, we have found that some chan- 
nels exhibit resonance states, in which the spin density 
contribution at the Fermi energy is extremely large (up 
to ~ 50 times larger than "typical" channels). In order 
to properly account for such states' contributions to the 
spin density, we have isolated such points in the Brillouin 
zone (if their contribution to the spin torque/current is 
greater than 2 fis/e - such states constitute typically 
at most 1% of the Brillouin zone area), and integrated 
over energy from to .8 mV (the potential required to 
reach typical critical current densities). This smooths 
out their contribution to the non-equilibrium spin den- 
sity. We have verified that points in the Brillouin zone 
which do not exhibit such resonances at the Fermi energy 
contribute an essentially constant spin density over this 
energy range, and may be safely treated within linear 
response. 

Once a self-consistent coUinear calculation is com- 
pleted, we initialized the various non-coUinear systems 



by rotating the initial spin density of the free layer to 
the desired angle. We found that non-coUinear config- 
urations do converge to self-consistency, and ostensibly 
represent extrema of the total energy. This is perhaps 
surprising since the exchange coupling between the two 
magnetic layers implies that the energy is extremized 
in parallel or anti-parallel configurations only. We have 
calculated the difference in energy between parallel and 
anti-parallel alignments and have found the difference to 
be less that 10^^ Hartree/orbital, thus verifying that to 
within the tolerance we use for self-consistency (a maxi- 
mum change of 10~^ Hartree in the Hamiltonian between 
iterations of the self-consistent cycle), a noncoUinear con- 
figuration can be a self-consistent solution. Our value for 
the change in energy between parallel and anti-parallel 
(or the exchange energy) is consistent with previously 
found values 



V. RESULTS 

Fig. 1 shows the calculated conductance per area ver- 
sus angle. From the conductivity in the parallel and anti- 
parallel alignment, a GMR ratio of 48% is obtained. The 
absolute value of the conductance is consistent with pre- 
vious calculations of similar systemsJ^ii^ The conduc- 
tance a.t 9 = 90° is not the mean between the conduc- 
tance of the parallel and anti-parallel state, which indi- 
cates that the behavior can not be captured by a simple 
rotation of the coUinear transmission coefficients in spin- 
space, but also includes higher order effects such as mul- 
tiple spin-dependent scattering between the two layers. 

Fig. 2 shows the spin torques as a function of the lay- 
ers' relative orientation. An important distinction is that 
between in-plane torques, which are non-energy conserv- 
ing and present only in non-equilibrium cases, and out- 



5 



of-plane torques, which are the result of itinerant electron 
exchange, and are responsible for RKKY-like interactions 
between layers. (In the following we refer to the out-of- 
plane torques as STT(6').) According to Fig. 2, the ratio 
of these two torques varies with angle, but we find it to be 
consistently above 10%, which can have important impli- 
cations for the behavior and stability of the magnetiza- 
tion dynamics. Interestingly, we also find that the out- 
of-plane torque undergoes a sign change near 6 = 180°. 
Previous studies have found a difference in sign between 
these torques under certain conditions but the relative 
sign is usually constant over all angles. The presence 
of an angle dependent relative sign of the torques may 
also have interesting consequences for the dynamics of 
the magnetization. 



Spin torque per current vs. layers' relative orientation 
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FIG. 2: Torque per current versus relative layer orientation. 

The angular dependence of the conductance is partially 
responsible for the departure from perfect sine behavior 
in STT(0). This departure is encapsulated in g{9), de- 
fined as STT(6')/(I • sin(0)), which has been calculated 
for simpler models,'^ and is an extremely important pa- 
rameter in spin torque physics. It represents the amount 
of spin torque delivered per electron in the current flow, 
and is therefore a measure of the efficiency of spin trans- 
fer torque. Fig. 3 shows our calculated g{9), and, for 
comparison, that found in the original calulation of Slon- 
cewski. Note that the y-axis for the two curves are differ- 
ent by a factor of 2, so that at small angle, our calculated 
result is more than a factor of 2 smaller than that of Slon- 
czewski. That is to be expected, as Slonczewski's model 
calculation considered the limiting case of pure spin fil- 
ters, where minority spins are completely reflected, and 
majority spin completely transmitted through the mag- 
netic interfaces. In this sense the Slonczewski can be 
considered to give the most optimistic spin transfer effi- 
ciency. Interestingly our result shows a smaller difference 
in efficiency between 9 — and 9 = n as compared to 



g(9) for various models 
-Slonczewski 




0.75 



LU 



0.5 I 



^0.25 



30 



60 



90 
Angle 



120 



150 



FIG. 3: Angular prefactor of spin transfer torque g{6). Note 
the Slonczewski scale differs from our cafculated g{9) scale by 
a factor of 2. 



Slonczewski, which is more consistent with experimentai 
data.30 

One can make further comparison to experiment in 
finding the amount of torque delivered per current. This 
quantity can be extracted from experimental data by 
finding the slope of the linear relation of critical current 
versus applied field. The experiment in Ref. |23 gives a 
slope of .29 mA/.l T. The resulting spin torque per cur- 
rent (or spin torque efficiency g(9) at 9 — 0) is given 

by "^^2^A^°^ ~ Sbfj-B/e. Here a is the bulk mag- 
netic damping of Co, assumed to be .007, Mg is the bulk 
magnetization of Co, and V is the volume of the free 
layer. The value obtained for the efficiency in this way 
is nearly identical to that determined in point contact 
experiments.^^ This contrasts with our calculated g(9) of 
.11. Our calculated efficiency is smaller than experimen- 
tal by about a factor of 3. Remarkably, even Slonccwski's 
efficiency is smaller than that seen in experiment by 50% 
(assuming polarization of Co to be .4). Some diffusive 
modeis of spin transfer typically predict much smaller 
efficiencies, with values for Co-Cu structures of about 
1%.'^^ At first glance then, it is a puzzle as to why the effi- 
ciencies seen experimentally are so high. One potentially 
important consideration is that the deduced experimen- 
tal efficiencies are not directly measured, but are rather 
model dependent on the mechanism for switching (usu- 
ally simple coherent switching is assumed), and depend 
on parameters such as the magnetic damping, whose val- 
ues and details are not well known. 

Fig. 4 shows the layer resolved spin torques. As 
expected, the torques show a generally oscillatory de- 
cay, which is understood as the result of averaging over 
many transverse channels' oscillatory contributions.— 
The length of the free layer is 2.67 nm, thus our result 
that the transverse spin density is close to complete decay 
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is consistent with other approaches, in which the decay 
length of transverse spin in Co is found to be 3nm.— 

We estimate the propensity (or lack thereof) for a spin 
wave or noncoUinear structure in the free layer due to 
the non-uniform torques by comparing the magnitude of 
these torques (or effective fields) and the exchange field 
an atom experiences from its neighbors. The effective 
Heisenberg nearest neighbor coupling constant for bulk 
Co has been computed to be 1.085 mRyi^ The result- 
ing exchange constant between planes (considering only 
nearest neighbor interactions) is then 4.34 mRy. This 
corresponds to a field of about lO^T. For a current of 
1 mA, the figure shows that typical non-equilibrium ex- 
change fields to be on the order of .0075 T. The resulting 
deflections of the spin are much less than 1°. This slight 
departure from coUinearity indicates that the surface- 
effect aspect of spin transfer torques (more generally, 
non-uniform torques) are not manifest in thin layers. 



Relative orbital contribution to STT 
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FIG. 4: Breakdown of relative orbital contributions to spin 
torque in the ^-direction, for the 9 = 90° case. 

Fig. 5 shows the relative contributions to the spin 
transfer from the s, p, and d orbitals. Such a partition ig- 
nores hybridization between orbitals, but one can never- 
theless extract meaningful results from such a division.— 
We find the d-orbitals' contribution is dominant. This is 
due to the substantial amount of d-electron conduction 
in the system, and more importantly to the relatively 
larger spin-dependent exchange-correlation potential of 
the d-electrons. A commonly used model for spin torque 
calculations in multilayers or in domain walls is the s-d 
model, in which the d-electrons are responsible for the ex- 



change field, while the s-electrons are responsible for the 
current. Spin torque on the d-electrons' magnetization is 
a result of the interaction of these two subsystems. How- 
ever, in the ballistic limit for Sd transition metals, our 
results show that it is largely the interaction between 
equilibrium and non-equilibrium d-electrons that is re- 
sponsible for spin torque physics. 



VI. SUMMARY AND CONCLUSIONS 

In this work, we have formulated a general scheme 
for calculating the dynamics of an order parameter in 
the presence non-equilibrium, current carrying quasipar- 
ticles. This scheme is applied to find the current-induced 
torques that are present in magnetic spin valve structures 
under bias. For the specific system studied, we have 
found the STT to be localized in the surface planes of 
the free magnetic layer, and have found the out-of-plane 
torque to be a substantial fraction of the in-plane torque. 
The magnitude of the torque is appreciably smaller than 
that deduced from experiment, although we note that 
experimental measures of the absolute spin torque are 
model-dependent. In addition, we have found the STT to 
be due mostly to the interaction between non-equilibrium 
d-electrons with the exchange field, which is itself also 
due mostly to the spin splitting of the d-band electrons. 

The formalism presented here is more general than that 
typically used in calculations of current induced torques 
in magnetic materials. In particular, since it does not 
rely on conservation of angular momentum, it may be 
applied to system with spin-orbit coupling, or to antifer- 
romagnetic systems. And the calculation of current in- 
duced torques on an atomic, or even atomic orbital scale 
resolution allows for the study these effects in molecu- 
lar scale systems. More generally, it may be applied to 
any mean-field system with an order parameter - for ex- 
ample superconductors, or ferromagnet-superconductor 
hybrids. 
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(AxS)/l - spin current-induced exchange torque and magnetic field per current 
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FIG. 5: Calculated layer resolved torque per current density and spin transfer induced change field per A. The (X,Y,Z) directions 
in the legend refer to the spin torques, while the direction of the effec 
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